Predicting and verifying transition strengths from weakly bound molecules 
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We investigated transition strengths from ultracold weakly bound 41 K 87 Rb molecules produced 
via the photoassociation of laser-cooled atoms. An accurate potential energy curve of the excited 
state (3) 1 E + was constructed by carrying out direct potential fit analysis of rotational spectra 
obtained via depletion spectroscopy. Vibrational energies and rotational constants extracted from 
the depletion spectra of v' =41-50 levels were combined with the results of the previous spectroscopic 
study, and they were used for modifying an ab initio potential. An accuracy of 0.14% in vibrational 
level spacing and 0.3% in rotational constants was sufficient to predict the large observed variation 
in transition strengths among the vibrational levels. Our results show that transition strengths from 
weakly bound molecules are a good measure of the accuracy of an excited state potential. 

PACS numbers: 34.20.-b,33.15.Bh,37.10.Mn,82.80.Ms 



I. INTRODUCTION 

Ultracold molecular gas is a prominent candidate for 
realizing novel distinctive applications in physics and 
chemistry, including precision measurements, quantum 
computation, ultracold chemistry, and novel quantum 
phases 0, Thus far the production of ultracold 

molecules in the vibrational ground state has been depen- 
dent on the optical transfer of weakly bound molecules 
formed via either photoassociation 0] or magnetoassocia- 
tion [i[ • Among the previously proposed methods IB-flol. 
stimulated Raman adiabatic passage (STIRAP) (ll| - Fl3j 
of weakly bound molecules is the most efficient method 
for preparing a molecular sample in a single quantum 
state. 

In the STIRAP transfer, an excited state is used as 
an intermediate state. It is important to select an ex- 
cited state having large transition strengths from both 
weakly bound and deeply bound levels. In general, the 
strengths of molecular transitions from low vibrational 
levels are readily predicted on the basis of the Franck- 
Condon factors (FCFs) calculated from a potential en- 
ergy curve or molecular constants 14] . The small number 
of nodes in the radial wavefunction for a low vibrational 
level indicates a low sensitivity to deviations in the ra- 
dial direction, thereby enabling us to explain the inten- 
sity distribution over vibrational levels from molecular 
constants. However, it is difficult to predict transition 
strengths from weakly bound molecular levels because 
weakly bound levels have a large number of nodes in their 
wavefunctions; hence, the FCFs are quite sensitive to the 
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FIG. 1. Relevant potential energy curves of KRb. Weakly 
bound molecules were produced by the photoassociation of 
laser-cooled 41 K and 87 Rb atoms. The (S) 1 ^ 4 " state was in- 
vestigated via the spectroscopy of molecules in the v" = 91 
level of X S + . The intermediate state for REMPI ionization 
was the (4) 1 E + state. 



wavefunctions of excited states. In this letter, we show 
that transition strengths from weakly bound levels can 
be predicted on the basis of an accurate potential energy 
curve of the excited state, constructed via direct poten- 
tial fit (DPF) analysis [l5[ of both vibrational energies 
and rotational constants. 

We focused on the (3) 1 S + state of KRb in the range 



2 



11400-11800 cm -1 with respect to the ground atomic 
threshold, which was proposed as a potential candidate 
for the STIRAP transfer of weakly bound molecules to 
the rovibrational ground state An RKR potential 

curve of the (3) 1 S + state was reported from the potential 
minimum up to 10400 cm" 1 [17| . An ionization spectrum 
for 39 K 87 Rb obtained by using a pulse laser and a deple- 
tion spectrum near 11700 cm -1 is provided in Ref. Ilq , 
Recently, we have realized the STIRAP transfer of weakly 
bound molecules («" = 91, J" = of A 1 E+) to the 
rovibrational ground state (v" — 0, J" — of X 1 E + ), 
mediated by the v' = 41 level of the (S) 1 !^ state [i~8j| . 
Before conducting this experiment, we carried out deple- 
tion spectroscopy in the range 11400-11800 cm" 1 , which 
revealed 10 vibrational levels. We found that the width 
of the observed spectra was highly dependent on the vi- 
brational levels of the (3) 1 S + state. By introducing an 
analytical representation for power broadening, we ex- 
tracted the transition strengths for each vibrational level. 
The other electronic states that correlate with the (3) 1 E + 
state via spin-orbit interaction were far away; hence we 
assumed that most of the perturbations from these elec- 
tronic states were negligible. Thus, we could analyze the 
experimentally obtained spectra on the basis of a sin- 
gle potential curve. In addition to the previous spectro- 
scopic work near the bottom of the potential (l7| , vibra- 
tional energies and rotational constants extracted from 
the spectra were used to construct an accurate potential 
via DPF analysis. Using the modified potential, it was 
possible to explain the variation in transition strengths 
among the vibrational levels of the (3) 1 E + state in terms 
of the FCFs. 



II. DEPLETION SPECTROSCOPY 

Previously, our experimental setup for the spec- 
troscopy of ultracold 41 K 87 Rb molecules was described 
in detail [19j. We provide a brief summary herein. We 
started with a dual-species magneto-optical trap (MOT) 
of 1 x 10 8 41 K atoms and 2 x 10 s 87 Rb atoms. A com- 
pressed MOT (C-MOT) procedure was applied for 40 ms 
to compress and cool the 41 K and 87 Rb atoms. The typi- 
cal densities and temperatures of 41 K and 87 Rb at the end 
of C-MOT were 2 x 10 11 cm" 3 and 400 /jK and 4 x 10 11 
cm" 3 and 100 /iK for 87 Rb, respectively. A photoassoci- 
ation (PA) laser (wavenumber, 12570.13 cm -1 ; intensity, 
1 x 10 3 Wcm" 2 ) was applied for 10 ms at the end of 
the C-MOT process. The produced molecules were de- 
tected using micro-channel plates (MCP) after they were 
ionized via resonance enhanced multi-photon ionization 
(REMPI) using a pulsed dye laser (wavenumber, 16543 
cm" 1 ; intensity, 3xl0 6 W cm" 2 ). 

Depletion spectra were obtained by monitoring ion 
counts in the v" = 91 level of A 1 S + , while a 
CW TkSapphire laser (Sirah Matisse TX; intensity, 50 
Wcm -2 ; beam waist, 350 fim) was continuously applied 
and scanned. In the present study, we analyzed spec- 
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FIG. 2. Depletion spectrum of v' = 41 level of (S)^" 1 taken 
with molecules in v" = 91 level of X X E+. The spectrum 
gives rotational constants of both ground and excited states. 
In addition, the width of the spectrum enables us to estimate 
the transition strength. 



tra for the v" = 91, J" = 2 level, whose binding en- 
ergy with respect to the atomic threshold FK=l+FRb=l 
was measured as —12.454(1) cm" 1 . The frequency of 
the TkSapphire laser was monitored using a Fabry-Perot 
cavity which was locked to 87 Rb D2 line. The cavity 
transmission signal was used to calibrate the variation in 
the scanning speed and to measure the relative frequency 
of the laser with respect to the cavity transmission peak 
with a precision of 3 MHz. The absolute frequency was 
measured using a commercial wavemeter (Wavelength 
WS-7; accuracy, 60 MHz). The accuracy of measure- 
ments for rotational constants was limited by the spectral 
width, which was of the order of 100 MHz, whereas that 
for vibrational energies was limited by both the spectral 
width and the wavemeter. 

Fig- 12 shows a depletion spectrum for the v' — 41 level 
of the (3) 1 S + state, which was used as an intermedi- 
ate state for the STIRAP transfer from the v" = 91, 
J" = level of A 1 S+ to the v" = 0, J" = level [H|. 
We can extract transition strengths as well as rotational 
constants for both ground and excited states from the 
spectra. In the appendix, we show that an approximate 
representation of the full-width-half-maximum (FWHM) 
of a depletion spectrum is given by 



2A - 0.79f2Vrr. 



(1) 



where is the Rabi frequency; T, the natural width of 
the excited state; and r, the duration of spectroscopy. 
Roughly speaking, the width increases not only with the 
light intensity but also with the duration of spectroscopy. 
In our case, the duration was estimated as a few millisec- 
onds on the basis of the temperature of the molecules and 
beam diameter of the depletion laser. The natural width 
of the (3) 1 S + state was not precisely known, but it was 
obtained as 27r x 300 kHz from an ab initio calculation 
[20| . In the following discussion, we derive the Rabi fre- 
quency from the observed spectrum by assuming r as 2 
ms, r as 2-7T x 300 kHz, and vTV as 60. 
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FIG. 3. Comparison between calculated and observed ro- 
tational constants. Calculated rotational constants of five 
molecular states in this range are plotted against energy levels 
from the ground atomic threshold. The observed rotational 
constants are in good agreement with those obtained from 
the ab initio potential of the (3) 1 E + state. There is a sys- 
tematic deviation of ~1% between the observed values and ab 
initio values, which indicates the inaccuracy of the ab initio 
potential. 
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FIG. 4. Potential energy curves of state. The modified 

potential obtained in the present study is compared with the 
ab initio potential and the experimental RKR potential. The 
shaded area denotes the region where depletion spectra were 
obtained. 



III. ANALYSIS 

Although the ab initio potential [2l[ enabled us to 
identify the symmetry of the observed spectra without 
any ambiguity, it gave rotational constants that exceeded 
the experimental ones by approximately 1% (Fig. [3]). 
There were two reasons for this deviation. First, the 
potential minimum exceeded the minimum of the RKR 
curve by 97 cm -1 . Second, the outer turning point was 
at the shorter internuclear separation than the RKR po- 
tential. We found that the variation in the observed tran- 
sition strengths could not be understood with the ab ini- 
tio potential. The observations could not be attributed 
to the small change in the vibrational quantum number, 
which corresponded to the energy difference in the min- 
imum of the RKR and the ab initio potentials. These 
facts indicate that the ab initio potential is inaccurate. 
In order to obtain an accurate potential, we carried out 
DPF analysis, whereby a potential is iteratively modi- 
fied until its eigenvalues coincide with those determined 
from the experimental spectra [HJ]. With the aid of the 
phiFIT program code |22|, we first constructed an an- 
alytical Extended Morse Oscillator (EMO) potential of 
the (3) 1 E + state on the basis of (1) the RKR curve, (2) 
a few points from the inner curve of the ab initio po- 
tential, and (3) a few points around 12500 cm -1 from 
the outer curve of the ab initio potential. The potential 
curve of the ground state A 1 E + was also required for the 
calculations. We used an EMO potential fitted to an ac- 
curate, experimentally determined potential [23| . Then, 
the analytical potential of the (3) 1 E + state was modified 
to reproduce our data by using the DPotFit program 
code [24j. A good convergence was achieved when we 



modified the potential significantly by manually moving 
the points from the ab initio potential. This procedure 
was repeated until the eigenvalues of the potential were 
within 0.05 cm -1 of the observed levels, i.e., only 0.14% 
of the vibrational level spacing. The remaining devia- 
tions were presumably due to the incomplete analytical 
function used to represent the potential. Table Q] lists the 
final potential parameters. These parameters are used 
to represent the potential energy curve in the following 
form: 

V(R) =K nin + D e (l- e -^)(«-«=)) 2 
12 

i=0 

R 3 



y(R, R 



Rt 



R 3 + R 3 ' 



(2) 



Fig. [4] shows the RKR potential, ab initio potential 
and modified potential. At the new potential, rotational 



TABLE I. Parameters for EMO potential obtained via DPF 
analysis of the observed spectra. The units are cm -1 for V m in 
and D e , Afor R e , and A -1 for 0j. 



Parameter 


Value 


Parameter 


Value 


Vmin 


9777.6963 


05 


-0.804673587 


D e 


3246.0363 


06 


3.00474008 


Re 


5.25904119 


07 


11.3877324 


00 


0.449794066 


08 


-2.140717462 


0i 


0.200265883 


09 


-39.73110586 


02 


0.406840126 


010 


-17.01400311 


03 


0.207884795 


011 


50.93279647 


04 


-0.476349301 


012 


44.83137419 
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TABLE II. Comparison between observed, ab initio, and modified values for vibrational energies (E) and rotational constants 
(B). Vibrational energies are energy levels of the J' = 1 level. Error bars in vibrational energies are 2 x 10 -3 cm" 1 , whereas 
those in rotational constants are 5 x 10 -5 cm" 1 . Observed vibrational energies are reproduced within 5 x 10 -2 cm -1 , whereas 
observed rotational constants are reproduced within 7 x 10 -5 cm -1 . The vibrational numbering for ab initio values is deviated 
by 2 because the potential minimum lies below that of the correct potential. 
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FIG. 5. Observed and calculated transition strengths. Ob- 
served values are determined from the spectra obtained via a 
depletion laser at 100 mW (red circle) and 500 mW (blue rect- 
angle). Calculations based on the modified potential can ex- 
plain the variation in transition strengths with respect to the 
vibrational levels in terms of the Franck-Condon factors. The 
point deviated from other points corresponds to the v' — 44 
level, which is expected to be mixed with other states. The 
inset shows the same plot based on the ab initio potential. 



constants were reproduced within 0.3% of the observed 
values (Table I3J. We excluded the v' = 44 level from 
the analysis because the characteristics of this level are 
anomalous; a much larger transition strength, a much 
smaller rotational constant, and a much larger deviation 
in a vibrational energy were observed at this level. These 
features indicate that this level was coupled to the (2) 1 E + 
state which was observed at ~1 cm -1 above the v' = 44 
level. 

The most important implication of this analysis is that 
the variation in the transition strengths with respect to 
the vibrational levels of the excited state can be accu- 



rately explained on the basis of the new potential. Fig. [3] 
shows a plot of the Rabi frequencies derived from the 
observed spectra against those calculated from the cor- 
rected potential. For comparison, the same plots based 
on the ab initio potential are also shown. The calcu- 
lated values are calibrated on the basis of our recent 
measurement of the transition dipole moment between 
(3) 1 E+, v' = 41 and X 1 ^ , v" = 91, 0.035(2)ea , de- 
termined via dark resonance spectroscopy of the rovibra- 
tional ground-state molecules. Weakly bound molecules 
have more than 90 nodes in their wavefunction; hence, 
the FCFs are highly dependent on the wavefunction of 
the excited state. In other words, the FCFs can serve 
as a sensitive measure of the accuracy of wavefunctions. 
Our results show that accuracies of 0.14% in vibrational 
level spacings and 0.3% in rotational constants are suf- 
ficient to predict the FCFs from weakly bound levels; 
these values arc justified by considering the typical size 
of nodes in the radial wavefunction. On one hand, the 
weakly bound level v" = 91 in the ground state X 1 ^ 
has an outer turning point of ^10 Aand an inner turning 
point of ~2 Ain the internuclear distance. Within these 
two points, there are 91 nodes; hence each node has a typ- 
ical size of ^0.1 A. Therefore, the required accuracy for 
representing the wavefunction is ~10 -2 A. On the other 
hand, in the present analysis, an accuracy of ^0.3% in 
rotational constants or ^0.15% in internuclear distance 
is obtained for the (3) 1 S + state because the relation be- 
tween the rotational constant B and the internuclear dis- 
tance R is given by B cx i?" 2 . Assuming the typical size 
of molecules in the v' =41-50 levels of (3) 1 S + as 6 A, 
we can derive the accuracy of the modified potential in 
the radial direction as ^10" 2 A; this value is in good 
agreement with the required accuracy for representing 
the weakly bound level. 

Now that we obtained an accurate potential curve as 
well as the absolute values of the transition dipole mo- 
ment of the v' = 41, (3) 1 E + level with the v" = and 
v" = 91 levels of we can predict the transition 
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FIG. 6. Transition dipole moments of the (3) E + state with 
the X E + state calculated from the experimentally obtained 
accurate potential. The absolute values for weakly bound 
levels are calibrated by dark resonance spectroscopy for the 
transition v' = 41, (3) X E + «- v" = 91,X 1 E+ whereas those 
for the rovibrational ground state are calibrated by dark reso- 
nance spectroscopy for the transition v' = 41, (3) 1 E + <— v" = 
0,X 1 E+. 



FIG. 7. Energy levels relevant to depletion spectroscopy 
of photoassociated molecules. Weakly bound molecules are 
formed at a rate of A pa- A laser having Rabi frequency 
of detuned from the resonance by A excites ground state 
molecules. Excited molecules spontaneously decay to the ini- 
tial level and other levels at rates of Ti and F2 — Ti, re- 
spectively. We mainly consider a situation where T2 > Vi is 
satisfied. 



dipole moment for each transition. Fig. [6] shows our pre- 
diction for the transition dipole moments of the (3) 1 S + 
state with the least bound state and for those with the 
lowest rovibrational level (v" = 0). The v' — 41 level 
used in Ref. [IH has favorable wavefunction overlaps with 
both the weakly bound and the lowest rovibrational lev- 
els; however other levels such as v' = 38 and v' = 39 can 
potentially serve as an intermediate state for the STIRAP 
transfer of weakly bound molecules to the rovibrational 
ground state. The potential presented herein can enable 
an accurate prediction for other isotopic combinations 
of KRb. Further, the present method for achieving an 
accurate potential and verifying its accuracy can be ex- 
tended to other molecular states that exhibit significant 
spin-orbit mixing by evaluating eigenvalues via coupled 
channel calculations including spin-orbit interaction. 



pare the transition strengths for each vibrational level. 
Rotational constants extracted from the observed spec- 
tra showed a 1% deviation from those calculated using 
an ab initio potential. By carrying out DPF analysis, we 
constructed an accurate potential that reproduced en- 
ergy levels with an accuracy of 0.14% in vibrational level 
spacing and 0.3% in rotational constants. The variation 
in transition strengths among vibrational levels could be 
understood in terms of FCFs calculated with the mod- 
ified potential. Our results indicate that the transition 
strengths from weakly bound levels serve as a sensitive 
measure of wavefunctions, which can be used to test the 
accuracy of the potential curve. In general, the proposed 
procedure can be adopted for constructing an accurate 
potential and verifying its accuracy on the basis of rota- 
tional spectra for weakly bound molecules. 



IV. CONCLUSION 



The (3) 1 S + state of KRb was investigated via the de- 
pletion spectroscopy of ultracold molecules formed by 
the photoassociation of laser-cooled 41 K and 87 Rb atoms. 
The spin-orbit mixing of other electronic states with the 
(3) 1 S + state was negligible; hence, we could assume this 
state as a single potential. The simplicity of the (3) 1 E + 
state enabled us to modify the potential to reproduce 
our observations as well as to assign the spectra. We 
observed 10 vibrational levels in the range 11400-11800 
cm -1 with respect to the ground atomic threshold. We 
developed a theoretical model that related the spectral 
width with the Rabi frequency, which was used to corn- 



Appendix: Analytical expression for the line shape 
of a depletion spectrum 

We consider a general situation, as shown in Fig. [3 and 
we assume that radiative transitions occur much faster 
than the time variation of the population in each molec- 
ular level because of the photoassociative creation of 
molecules. We first consider the evolution of the popula- 
tion in the ground and excited states of a single molecule. 
The optical Bloch equation for this system can be written 
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ds 

~dt 

dw 

~dt 
dt 



=-(r 2 -r x 
=(r 2 + r x ) 



w 



2 

s — w 



2illm(p eg ) 



T 2 i 

-(— -iA)p eg + -wCl 



(A.l) 



where s — p gg + p ee , w = p gg - p ee , and |e > and \g > 
denote an excited state and a weakly bound level, re- 
spectively. A and 51 denote the detuning frequency and 
the Rabi frequency, respectively. The decay rates from 
the excited state to the initial weakly bound level and 
to other levels are given by I\ and T 2 , respectively. The 
photoassociation rate is denoted by Ap^. The width is 
much larger than the Rabi frequency in the experiment; 
hence, we can assume that the time evolution of s and 
w is much slower than that of pT g . Thus, we can set 
dp^ g /dt — and obtain the following expression for p^ g : 



Peg 



T 2 - 2iA 



(A.2) 



Substituting Eq. (|A~2l) in Eqs. (lA~Tj) . we obtain alternative 
equations for s and w as 



dw -rr + r ) s ~ w 



T 2 + 4A 2 



(A.3) 



Here, we assume that the time taken by the mean value of 
the population ratio w/s to attain a constant value z after 
is greater than the typical time for radiative transitions. 
By using the relation w = zs, Eqs. (|A.3p give following 
equations: 



ds , > 1 - z 

ds ,_, - .1 — z 



rl + 4A 2; 



(A.4) 



Substituting the first equation in the second equation, we 
obtain a time-independent equation for z: 



(1-z) \T 3 (l + z)+T 1 (l-z) 



4» 2 r 2 
r 2 2 + 4A 2 ; 



(A.5) 



This equation is readily solved, and it gives the following 
expression for z. 



^/(l + k 2 )T 2 + 2T 1 T 2 k - (kT 2 + r t ) 



r 2 -r x 



k = 



2n 2 



r 2 + 4A 2 



(A, 



This expression is used in the following discussion. Next, 
we derive rate equations for the population in the ground 
and excited molecular levels for the number of molecules: 



dN 
~dt 



--A PA - (r 2 

l ' Z -N 



2 

1 + z 



N 



(A.7) 



where N g and N e denote the number of molecules in the 
ground and excited states, respectively, and = N g + N e 
is the total number of molecules. In these rate equa- 
tions, a typical timescale is of the order of 1 ms, and it is 
governed by A pa. This is much longer than the typical 
timescale for radiative transitions in most cases (<1 p,s). 
The time evolution of N g is given by 



dN n 



1 + z A 1 



dt 2 2 

Thus, the solution for N g is given by 

A PA l + Z 



- T\)N g 



(A- 



N g (r) = 



r 2 - ri i - 

I 1 — exp 



1 



: (r 2 -ri)T 



(A.9) 



This expression gives the line shape of a depletion spec- 
trum for a duration r. Assuming (T 2 — Ti)t ^> 1, 
Eq. (|A~9)) gives 



N g (r) -> 



A PA 

r 2 -r : 



(A.10) 



on resonance (A = and w — > 0), whereas N g at an 
infinite detuning (A — > oo and w — > 1) is given as 



N 9 (t) -> A PA T 



(All) 



N g at an infinite detuning is much larger than N g on res- 
onance. Thus, the width of the spectrum is determined 
by finding w such that it satisfies 



1 



N g (r) = -Ap A T 



(A.12) 



Substituting Eq. (|A.9|) in Eq. (IA.12p and rewriting the 
equation with a new variable x = 1 — z(<C 1), we ob- 
tain the following equation: 



1 — exp 



:(r 2 -r x )T 



(A.13) 



A rigorous solution of this equation is given by 
_4 + 2W(-2/e 2 ) 



(r 2 



i 



(r 2 -r!)r 



Fi)r 

x 3.18724... 



(A.14) 
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where W is the Lambert W function. Thus, the param- 
eter k in Eqs. (|A.6[) is given by 

_ 3.18724(1 - 1.59362/IV) 

(r 2 -r!)r- 3.18724 (A ' i5j 

Assuming (r 2 — Ti)r ^> 1, we obtain the following ex- 
pression for FWHM: 

2A«0.79fV(r 2 -ri)T ( A - 16 ) 

When decays from the excited state to the initial state 
are negligible (r 2 ^s> Ti), we obtain a simple relation 
([1]). It is difficult to evaluate the numerical factor yTr 
precisely; therefore, the width of a depletion spectrum 



cannot serve as an accurate measure of the transition 
strength. However, the expression (JTJ enables us to sys- 
tematically compare transition strengths for different vi- 
brational levels. 
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